Conformal invariance in two-dimensional turbulence 
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Simplicity of fundamental physical laws manifests itself in fundamental sym- 
metries. While systems with an infinity of strongly interacting degrees of free- 
dom (in particle physics and critical phenomena) are hard to describe, they often 
demonstrate symmetries, in particular scale invariance. In two dimensions (2d) 
locality often promotes scale invariance to a wider class of conformal transfor- 
mations which allow for nonuniform re-scaling. Conformal invciriance allows 
a thorough classification of universality classes of critical phenomena in 2d. Is 
there conformal invciriance in 2d turbulence, a peiradigmatic example of strongly- 
interacting non-equilibrium system? Here, using numerical experiment, we show 
that some features of 2d inverse turbulent cascade display conformal invariance. 
We observe that the statistics of vorticity clusters is remarkably close to that of 
critical percolation, one of the simplest universality cleisses of critical phenom- 
ena. These results represent a new step in the unification of 2d physics within 
the framework of conformal symmetry. 

We consider here 2d incompressible turbulent motion of a fluid, which represents an 
appropriate description of large-scale motions of the atmosphere and can be realized in 
different laboratory settings as well^~®. As predicted by Kraichnan^, stirring at some scale 
Lf results in two turbulence cascades, with the formation of fine-scale vortical structures and 
large-scale velocity structures. In two dimensions, squared vorticity uj"^ — (V x vY performs 
a direct cascade to small scales while kinetic energy flows from the injection length I// to 
large scales, opposite to the three-dimensional case. We focus here on the inverse cascade of 
energy for which, not surprisingly in view of the presence of a strong interaction, there is no 
exact analytic theory. Phenomenological dimensional arguments give consistent predictions, 
though in two seemingly unrelated ways. Consider the velocity difference at the distance 
r. On the one hand, one may require that the kinetic energy v'^ divided by the typical 
time r/vr must be constant and equal to the energy flux, e: ~ er. On the other hand, 
it can be argued that vorticity, which cascades to small scales, must be in equipartition in 
the inverse cascade range^. If this is the case, the enstrophy r'^uj^ accumulated in a volume 
of size r is proportional to the typical time r/vr at such scale, i.e. r'^uj'^ ~ r/v.f.. Using 
uOj. ~ v.,./r we derive ~ r^"*^ which for c? = 2 is exactly the requirement of constant energy 
flux. Amazingly, the requirements of vorticity equipartition (i.e. equilibrium) and energy 
flux (i.e. turbulence) give the same Kolmogorov-Kraichnan scaling in 2d. Experiments^'^'^ 
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and numerical simulations^ indeed demonstrate scale-invariant statistics with the vorticity 
having scaling dimension 2/3: cUr (x r~^/'^. 

Our goal here is to find out whether scale invariance can be extended to conformal 
invariance at least for some properties of 2d turbulence. Under conformal transformations 
the lengths are re-scaled non-uniformly yet the angles between vectors are left unchanged 
(a useful property in navigation cartography where it is often more important to aim in 
the right direction than to know the distance)^'^^. The novelty of our approach is that 
we analyze the inverse cascade by describing the large-scale statistics of the boundaries of 
vorticity clusters, i.e. large-scale zero- vorticity lines. In equilibrium critical phenomena, 
cluster boundaries in the continuous limit of vanishingly small lattice size were recently 
found to belong to a remarkable class of curves that can be mapped into Brownian walk 
(called Stochastic Loewner Evolution or SLE curves)^^"^^. Namely, consider a curve 7(t) 
that starts at a point on the boundary of the half-plane H (by conformal invariance any 
planar domain is equivalent to the upper half plane). One can map the half-plane H minus 
the curve 7(t) back onto H by an analytic function gt{z) which is unique upon imposing the 
condition gt{z) ~ z + 2t/z + 0{l/z'^) at infinity. The growing tip of the curve is mapped into 
a real point C,{t). Loewner^*^ found in 1923 that the conformal map gt{z) and the curve 7(t) 
are fully parametrized by the driving function ^(t). Almost eighty years later, Schramm^^ 
considered random curves in planar domains and showed that their statistics is conformal 
invariant if ^(t) is a Brownian walk, i.e. its increments are identically and independently 
distributed and ((^(t) — ^(0))^) = nt. In simple words, the locality in time of the Brownian 
walk translates into the local scale-invariance of SLE curves, i.e. conformal invariance. SLE^ 
provide a natural classification (by the value of the diffusivity k) of boundaries of clusters of 
2d critical phenomena^^ described by conformal field theories (CFT)^'^ and allow to establish 
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many new results (see 

The fractal dimension of SLE^ curves is known to be^^'^^ = 1 + k/8 for /t < 8. To 
establish possible link, let us try to relate the Kolmogorov-Kraichnan phenomenology to the 
fractal dimension of the boundaries of vorticity clusters. Note that one ought to distinguish 
between the dimensionality 2 of the full vorticity level set (which is space-filling) and a 
single zero- vorticity line that encloses a large-scale cluster^^. Consider the vorticity cluster of 
gyration radius L which has the "outer boundary" of perimeter P (that boundary is the part 
of the zero-vorticity line accessible from outside, see Fig. Olfor an illustration). The vorticity 
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flux through the cluster, / udS ~ ujlL"^, must be equal to the velocity circulation along the 
boundary, F = / v-d£. The Kolmogorov-Kraichnan scaling is ujl ~ (coarse-grained 
vorticity decreases with scale because contributions with opposite signs partially cancel) so 
that the flux is oc L^/'^. As for circulation, since the boundary turns every time it meets 
a vortex, such a contour is irregular on scales larger than the pumping scale. Therefore, 
only the velocity at the pumping scale Lf is expected to contribute to the circulation, such 
velocity can be estimated as {eLfY^^ and it is independent of L. Hence, circulation should 
be proportional to the perimeter, F oc P, which gives P oc L^/^, i.e. the fractal dimension 
of the exterior of the vorticity cluster is expected to be 4/3. 

Let us check this hypothesis by data analysis. A powerful tool for the study of 2d 
turbulence is the numerical integration of the incompressible Navier-Stokes equations in a 
planar domain. By this method it is possible to achieve a range of dynamical lengthscales 
of about four decades, whereas current laboratory experiments are limited to a scale sepa- 
ration of about a hundred. We present here the analysis of very high resolution numerical 
simulations (with up to 16384^ grid points) of two-dimensional inverse cascade. (See Ta- 
ble H] for the details). Vorticity clusters are shown on Figure E The fractal dimension of 
their exterior boundary (without self-intersections) shown in the left panel of Figure El as 
indeed close to = 4/3. Moreover, the fractal dimension of the boundary itself is close 
to D = 7/4. Of course, having some particular dimension does not by itself imply that 
the curve belong to SLE. Note, however, that the exterior perimeter of SLE^ with k > 4 
is conjectured^^ to look locally as SLE^. curve with = 16/k < 4 resulting in the duality 
relation, (/} — 1)(D* — 1) = 1/4, as observed in our turbulence data. Moreover, =4/3 cor- 
responds to a SLEg/s curve which represents the continuum limit of a self-avoiding random 
walk while the dual SLEg curve corresponds to a cluster boundary in critical percolation. 
That prompts us to compare the probability distributions of sizes and boundary lengths be- 
tween vorticity and percolation clusters. The size s of a cluster is the number of connected 
sites with the same sign of vorticity, the boundary length b is the number of sites that belong 
to the cluster but are adjacent to sites of different sign, and the diameter L is the side of 
the smallest square that covers the cluster. The results shown in the right panel of Fig. |21 
are in a good agreement with the exact results from percolation theory. 

The two SLEs with k, = 6 and = 8/3 correspond to CFT with zero central charge 
which means that the scale invariance remains unbroken even when the system is on a 
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manifold with corners or with a nonzero Euler number (a topological invariant determined 
by the number of handles and boundaries). Also, SLEg curves are singled out by a "locality 
property" (the curve does not feel the boundary until it touches it), while their dual SLEg/s 
has the "restriction property" (the statistics of the curves conditioned not to visit some 
region is the same as in the domain without this region) How does all this relate to 

two-dimensional turbulence? 

Now we show by a straightforward check that, within statistical accuracy, large-scale 
zero-vorticity lines are indeed SLE curves, that is they are conformal invariant and possess 
remarkable properties of a kind that has never been studied in turbulence. Zero-vorticity 
isolines that are candidate SLE traces are identified as follows. First, a horizontal line 
representing the real axis in the complex plane is drawn across the vorticity field. Second, 
an explorer starting from the origin at the real axis walks on the zero-vorticity isoline keeping 
the positive vorticity sites always on the right. Third, when the explorer hits the real axis 
it treads on it always leaving the positive region on its right side until it can re-enter the 
upper half-plane. This eventually leads the explorer to infinity in an unbounded domain. 
An example of the outcome of this search is shown in Figure El Strictly speaking, this 
procedure faithfully reproduces all the details of the statistics only if there is a locality 
property (meaning that the exploration process does not feel the boundary before it hits 
it), which holds for SLEg. Since we obtain as a result SLE with k ^ 6, our procedure is 
self-consistent; we also checked that shifting and turning the line does not modify the results 
presented below. To determine which driving function ^(t) can generate such a curve, one 
needs to find the sequence of conformal maps gt{z) that map the half-plane H minus the 
curve into H itself. We approximate gt{z) by a composition of discrete, conformal slit maps 
that swallow one segment of the curve at a time (a slight variation of the techniques presented 
in This results in a sequence of "times" ti and driving values C,i that approximate the 

true driving functions. If the zero-vorticity isolines in the half-plane are actually SLE traces, 
then the driving function should behave as an effective diffusion process at sufficiently large 
times. We have collected 1,607 putative traces. The data presented by blue triangles in 
Figure m show that the ensemble average (^(t)2) indeed grows linearly in time: the diffusion 
coefficient k is very close to the value 6, with an accuracy of 5% (see inset). The average 
{C,{t)) vanishes by the reflection symmetry of the Navier-Stokes equations. Additionally, the 
probability distribution functions of ^ (t) / collapse onto a standard Gaussian distribution 
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at all times t. Therefore, we expect the driving ^{t) tend to a true Brownian motion and zero- 
vorticity lines to become SLE^ traces with k very close to 6 in the limit of vanishingly small 
Lf. To appreciate how remarkable this property is, pink symbols in the lower insert in FigE] 
show for comparison the results of the same procedure for the isolines of a Gaussian field 
having the same Fourier spectrum as vorticity but randomized phases. The slow incomplete 
recovery to k = 6 for the random-phase field occurs at the scales where the power-law 
correlation is already cut-off by friction and the field becomes truly uncorrelated. 

The identification of isovorticity lines as SLE^ curves allows to apply powerful techniques 
borrowed from the theory of stochastic differential equations and conformal mapping theory 
and to obtain analytic predictions for some nontrivial statistical properties of vorticity clus- 
ters. The first example is the probability that a point z = pe*^ inside the upper half plane is 
surrounded by a positive vorticity cluster connected to the positive real axis. In this event, 
it is not possible to reach infinity with a continuous path starting at z without treading 
on positive vorticity sites. For this to happen the zero- vorticity line must leave the point 
z on its right. The probability of such an event depends only on the angle 6 between the 
point and the origin and it assumes a particularly simple form in terms of hypergeometric 
functions^''. In the inset of Figure |S1 we show that the analytic solution fits very well the 
numerical data with n = 5.9. The second example is the probability that a vorticity cluster 
spans the rectangle joining two opposite sides. What is the average number of such spanning 
clusters? What is the probability that a "four-legged cluster" joins all four sides? By scale 
invariance these quantities depend only on the aspect ratio r of the rectangle, and their 
precise dependence can be found by exploiting conformal invariance. In the context of crit- 
ical percolation formulae for such probabilities have been derived by Cardy^^ and Watts^^ 
and later proven by Smirnov^^ and Dubedat^^. In the main frame of Figure El we show that 
numerical data for vorticity clusters follow very closely the expectations for SLEg. We have 
also checked (green symbols in Fig|2I) that the dimension of the set of narrow necks that 
enclose large fjords or large peninsulae has dimension^^ 3/4 (the set is defined by the pairs 
of points on the curve that are closer than Lf yet separated by an arclength larger than 
1000 Lj). All that gives further support to the result that zero- vorticity lines are conformal 
invariant and belong to the same class of universality as boundaries of percolation clusters. 

Whether the statistics of the zero-vorticity isolines indeed falls into the simplest uni- 
versality class of critical phenomena (and the fractal dimensions are exactly 7/4 and 4/3) 
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deserves to be a subject of more study. Do our findings signify that universal nature of 
percolation extends to turbulence as well as to diffusion-limited aggregation^^ and quantum 
chaos^^? At the present level it has the status of a tantalizing conjecture with strong - 
though not conclusive - support from the data. In view of the non-local constraint im- 
posed by the flow incompressibility, it is quite surprising that the statistics of zero-vorticity 
isolines (within experimental accuracy) enjoys the locality property inherited by its SLEg 
nature. Remind that continuous percolation can be constructed as a "flooded landscape" 
determined by some short-correlated random height function. However, vorticity field in 
the inverse cascade is not short-correlated, it has power-law correlation {uj(0)uj(r)) oc r~^^^. 
When the pair correlation function falls slower than r"'^/^ then the system is not expected 
generally to belong to the universality class of uncorrelated percolation and even be confor- 
mal invariant^^. Indeed, we have seen that the field having the same pair correlation function 
as the vorticity yet randomized phases of the Fourier harmonics does not have conformal 
invariant isolines (pink symbols in the lower inset in Fig|^. We thus conclude that there 
is indeed something special about the vorticity (which has nontrivial phase correlations and 
higher moments) produced by 2d turbulence. It is also intriguing to notice that conformal 
field theory of critical percolation possesses a field of scaling dimension 2/3, identical to 
the one for the vorticity in Kolmogorov-Kraichnan phenomenology. We may also wonder 
how conformal invariance is broken in statistical properties of non-zero vorticity isolines. 
Let us stress that we have found conformal invariance for zero-vorticity isolines, not yet for 
correlation functions as envisaged by Polyakov'^; in the related problem of passive scalar in 
turbulent flow correlation functions are not conformal invariant^'^^. 

To conclude, we have developed a numerical tool for testing conformal invariance in 
physical systems, established this symmetry (within experimental accuracy) for 2d inverse 
cascade and used it as a powerful new tool in turbulence study which allowed us to make new 
quantitative predictions confirmed by the experiment. That shows how conformal invariance 
spans the whole physics, from exalted subjects like string theory and quantum gravity, via 
statistical mechanics and condensed matter, down to earthly atmospheric turbulence. 
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FIG. 1: Vorticity clusters. These are defined as connected regions with the same sign of vorticity 
(here positive). Colours are arbitrarily attributed to different clusters. Regions of negative vorticity 
are black. The forcing lengthscale -^^/ is one hundredth of the box side. 



11 




B 

00 

o 



o 

CO 



10" 



10" 



^ 10" 



o 10" 



10" 



10" 




5 
4 
3 
2 
1 

10'^ 10""^ 10"^ 



1 


1 Q 






» 1 


1 



diameter L 



10"^ 10"^ 10""" 10"^ 
size (s) or boundary length (b) 



FIG. 2: Fractal dimensions and probabilities of size and boundary length for vorticity clusters. 
Left panel: The fractal dimensions are the slopes of the length-diameter dependencies in log-log 
coordinates for the boundary of filled clusters (red), for the outer boundary (blue) and for the 
necks of large fjords/peninsulae (green). The solid lines have slopes with the exact values for SLEg 
curves, 7/4, 4/3, 3/4, respectively. Fractal dimensions are obtained by computing the average 
length for a given diameter of the cluster. Right panel: The fraction of clusters with sizes between 
s and 1.25s (red symbols) and with boundary lengthes between b and 1.256. The solid lines are 
the predictions from the percolation theory. Inset shows the same data multiplied by s^^/gi g^.^^ 
6*/^, respectively. The vertical scale is linear to appreciate the plateau in the compensated plot. 
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FIG. 3: A portion of a candidate SLE trace obtained from the vorticity field. The red curve is a 
zero-vorticity hne in the upper half-plane. The dashed blue line is the "outer boundary" of the 
red curve, i.e. the boundary of the region that can be reached from infinity without getting closer 
than Lf to the red curve. The green dots mark the necks of large fjords and peninsulae. 
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FIG. 4: The driving function is an effective diffusion process with diffusion coefficent k = 6 it 0.3. 
The inverse cascade range corresponds to 5 • 10~^ < t < 10^^. Main frame: the hnear behaviour of 
(^(i)^). Lower-right inset : Diffusivity: blue for vorticity isoUnes, pink for the field with randomized 
phases. Upper-left inset: the probability density function of the rescaled driving function ^(t)/\/l^ 
at four different times t = 0.0012, 0.003, 0.006, 0.009; the solid line is the Gaussian distribution 
g{x) = (27r)-i/2exp(-xV2). 
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FIG. 5: Crossing and surrounding probability for vorticity clusters. Main frame: the probability vr^ 
that a cluster crosses from top to bottom a rectangle of aspect ratio r (blue), the average number A^^, 
of vertically crossing clusters (green), and the probability vr/^^ of a "four-legged" cluster joining all 
sides of the rectangle (red). The lines are the exact results for k = 6: vr^ = r]^'^ 2F1 ( 51 1; |; ^) 
with rj = [{1 - k)/{l + k)]^ and r = K{1 - k"^) / [2K (k^)] (Cardy-Smirnov, thick solid line); = 
\ [t^v + vr/i„ - log 7]] (Cardy, thick dashed line); -khv = t^v - Y{^)r{^) ^^"^ I' ^' I' ^) (Watts- 
Dubedat, thin dotted line). Inset: the probability that a zero-vorticity line in the upper half-plane 
leaves the point pe*^ to its right, for p = 0.048,0.064,0.080,0.096. The prediction for SLE^ traces 
is P = ^ -|- ^^j.^l.K.y 2F1 (^^) ^; |; — cot^ 9j cot 9, shown as a thick solid line for k = 5.9 (the best 
fit). The dashed lines are the probabilities for k = 5.7 and k = 6.1. 
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N dx V a Urms Lf 4 El Vi Vu Va 

2048 4.9 X 10"*^ 2 x 10"^ 0.015 0.26 0.01 2.4 x 10^^ 3.9 x 10"^ 1.8 x lO"''^ 2.1 x 10"^ 39.3 38.0 1.3 
4096 2.4 X 10-4 5 x 10"^ 0.024 0.26 0.01 1.2 x lO'^ 3.9 x lO'^ 0.7 x lO'^ 3.2 x lO'^ 39.3 36.1 3.2 
8192 1.2 X 10-4 2 X 10-6 0.025 0.27 0.01 7.8 x IQ-^ 3.9 x IQ-^ 0.3 x IQ-^ 3.6 x IQ-^ 39.3 35.3 4.0 
16384 0.6 X 10-4 1 X IQ-^ 0.0 0.24 0.01 5.5 x IQ-^ 3.8 x IQ-^ 0.2 x IQ-^ 3.6 x IQ-^ 39.5 37.6 1.9 

TABLE I: Parameters of the simulations. N spatial resolution, dx grid spacing, viscosity, a 
friction, Lf forcing scale, id = v^^'^/ijl^^ enstrophy dissipative scale, sj energy injection rate, Si, 
viscous energy dissipation rate, Sa energy dissipation by large-scale friction (energy growth rate 
for N = 16384), 77/ enstrophy injection rate, r)i, viscous enstrophy dissipation rate, rja enstrophy 
dissipation by friction (enstrophy growth rate for N = 16384). 
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